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ABSTRACT 

FN Tau is a rare example of very low-mass T Tauri stars that exhibits a spatially 
resolved nebulosity in near-infrared scattering light. To directly derive the parameters 
of a circumstellar disk around FN Tau, observations of dust continuum emission at 340 
GHz are carried out with the Submillimeter Array (SMA). A point-like dust continuum 
emission was detected with a synthesized beam of ~ 0.7" in FWHM. From the analysis 
of the visibility plot, the radius of the emission is estimated to be < 0.29", corresponding 
to 41 AU. This is much smaller than the radius of the nebulosity, 1.85" for its brighter 
part at 1.6 /im. The 340 GHz continuum emission observed with the SMA and the 
photometric data at A < 70 fj,m are explained by a power-law disk model whose outer 
radius and mass are 41 AU and (0.24 — 5.9) x 10~^ Mq, respectively, if the exponent of 
dust mass opacity (/3) is assumed to be — 2. The disk model cannot fully reproduce the 
flux density at 230 GHz obtained with the IRAM 30-meter telescope, suggesting that 
there is another extended "halo" component that is missed in the SMA observations. 
By requiring the halo not to be detected with the SMA, the lower limit to the size of 
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the halo is evaluated to be between 174 AU and 574 AU, depending on the assumed /3 
value. This size is comparable to the near-infrared nebulosity, implying that the halo 
unseen with the SMA corresponds to the origin of the near-infrared nebulosity. The 
halo can contain mass comparable to or at most 8 times greater than that of the inner 
power-law disk, but its surface density should be lower than that at the outer edge 
of the power-law disk by more than one order of magnitude. The physical nature of 
the halo is unclear, but it may be the periphery of a flared circumstellar disk that is 
not described well in terms of a power-law disk model, or a remnant of a protostellar 
envelope having flattened structure. 

Subject headings: circumstellar matter — submillimeter: planetary systems — proto- 
planetary disks — stars: individual (FN Tau) — stars: pre-main sequence 



Introduction 



Hundreds of exoplanets were discovered in the last 15 years (jMarcv et al.ll2005l : iMavor et al 



20051 ). Many planetary systems have a configuration completel y different froni the solar system, 
requiring us a comprehensive theory for planet formation (e.g., Ilda &: Linll2004l ). Diversity of the 
exoplanets may stem from a variety of the properties of their precursors, i.e., protoplanetary disks. 
Kokubo fc Idal (j2002l ). for example, proposed the idea that the planet configuration formed in 
a disk depends on its initial material distribution. Dust continuum emission at millimeter and 
submillimeter wavelengths is an important tracer to the total mass or surface density of the disks. 
Indeed, a sensitive multiwavel ength submillimeter c o ntinu um survey of 153 young stellar objects in 
the Taurus- Auriga regions by lAndrews &: Williamsl (|2005l ) revealed that circumstellar disk masses 
are lognormally distributed with a mean total mass of 5 x 10~^ Mq and a large dispersion (0.5 dex) 
when the gas-to-dust mass ratio of 100 is adopted. Meanwhile, imaging surveys of Sun-like single T 
Tauri stars were carried out with millimeter and submillimeter interferometers that can resolve disk 



emission (jKitamura et al.ll2002l : lAndrews &: Williamsl 120071 ). They obtained not only typical values 
or evolutionary trends of various disk parameters but also found scatters around them, which can 
be interpreted as a diversity. 

Typical values of disk parameters or what kind of planet is likely to form must also depend on 
the stellar mass. Planet formation around stars with a mass lower than Sun-like stars is especially 
intriguing, since these are the most abundant stars in our galaxy. IScholz et al.l (j2006l ) conducted 
sensitive 1.3 mm observations of 20 young brown dwarfs in the Taurus star- forming region and 
found that their disk masses range from < 0.4 to several Jupiter masses, indicating that the disk 
mass relative to the star is comparable to those derived for coeval low-mass (< 3 Mq) stars, 
namely, < 1 % to 5 % in most cases. They also found that at least 25% of the targets are likely 
to have disks with radii greate r than 10 AU, implying that the dynamical ejection of embryos 



(Bate. Bonnell. &: Bromm 



20021 ) might not be the dominant formation process for brown dwarfs 
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(see also lLuhman et al.l 120071 ) . These results may suggest that there may be no essential difference 
in the process for the formation of a star-disk system between Sun-like stars and less massive stars 
including brown dwarfs. Compared to the Sun-like stars, however, the typical mass and surface 
density of the disks around M dwarfs and brown dwarfs are expected to be lower while the snow line 
is located at a smaller radius. Such differences in the initial condition may s uppress the forrnation 



of gas giants bu t may result in higher frequency of icy Neptune-mass planets (jLaughlin et al 



2004 



Ida k. Linll2005l ). Althoug h these theoretical predictions seem to be supported by recent planet 



searches around M dwarfs (Beaulieau et al 



2006 



Johnson et al.l 120071 ). physical properties of the 



disks around their younger counterparts should be studied further to verify the above considerations. 

Another important clue for understanding the planet formation process in the disks is dust 
grain property such as composition, crystallinity, and size distribution. These have been studied by 
spectroscopic observations of silicate emission features at mid-infrared wavelengths and have been 
obs erved toward man y young stars ranging from brown dwarfs to Her b ig Ae /Be stars (for review. 



see 



Natta et al.l 120071 ). For Sun- like stars, for example, ISargent et al.l (120061 ) analyzed 8 — 14 fim 



emission spectra of 12 T Tauri stars in the Taurus- Auriga dark clouds and TW Hydrae obtained 
with Infrared Spectrograph (IRS) on board the Spitzer Space Telescope. They found that later 
spectral type stars can have relatively large amounts of crystalline silicates in their surrounding 
disks. More recently, a comparative study of the dust and gas properties of disks around young 



Sun- li ke stars (K1-M5) and cool stars including brown dwarfs (M5-M9) was made bv lPascucci et al 



([20091). They revealed that lower-mass stars including brown dwarfs tend to show weaker features 
at 10 ;um, indicating that the features are dominated by mor e "processed" du st grains than those 
from disks around higher- mass Herbig Ae/Be stars (see also lApai et al.ll2005l ). These systematic 
trends may be related to the fact that the location of the 10 //m silicate emission zone is closer 
to the star as the stellar luminosity is lower; the larger crystallinity seen toward cooler stars can 
be accounted for if the size of dust grains in the inner regions is systematically larger because the 
contribution to the spectrum from amorphous silicate s diminishes more rapidly with grain growth 
(0.1 — 10 /um) than that from crystalline silicates do (jKessler-Silacci et al.l 120071 ). Although these 



features arise from the optically thin surface layer of the inner regions of the disks and might not 
necessarily reflect the properties of the bulk silicate dust in the observed disk, these are the direct 
evidence for dust processing (i.e., dust growth and crystallization), which can be regarded as the 
first step to the formation of planets. 

FN Tau (Haro 6-2) is a relatively rare example of low-mass T Tauri stars whos e near-infrared 
scattering light shows clear appearance of disk- like morphology (jKudo et al.l 120081 ). It is a single 
star with a spectral type of M5, located at the distance of ^ 140 pc. Th e stellar mass and 



age are estimated to be 0.11 Mq and < 10^ yr, respectively ( MuzeroUe et al.ll2003l ). Existence 
of a circumstellar d isk was suggested by excess emission at infrared and millimeter wavelengths 



(| Strom et al.l 1 19891 : iBeckwith et al.l Il990l ) . Significant fraction of crystalline silicate grains was 



inferred from mid-i nfrared spectroscopic observations, as co mmonly be seen in the cases of other 



cool T Tauri stars (jSargent et al 



2006 



Forrest et al.l l2004l ) , suggesting that grain processing is 
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ta king place at least in the disk inner regions. Recent high-resolution imaging at A = 1.6 
by iKudo et al.l (j2008l ) discovered circumstellar nebulosity. The shape of the nebulosity is almost 
circular, and the azimuthally averaged radial profile of the brightness follows power-law forms of 
^-2.5±o.i HQ AU < r < 260 AU and r-6-5±0-2 280 AU < r < 400 AU. Based on the fact 
that the exponent of th e former law agrees with that observed in the nebulosity around TW Hya 
(jWeinberger et al.ll2002l). which can be explained i n terms of scattering light from the surface of 

they argued that the origin of near-infrared 



a flared disk (Krist et al 



200C 



Trilling et al 



2001 



nebulosity around FN Tau is also a flare d circumstellar disk. Assuming that dust emission at 1.3 
mm (31 lb 9 mJy; iBeckwith et al.lll99d ) comes from the same regions as the nebulosity, the disk 
mass was estimated to be 7 x 10~^ Mq. To derive the disk parameters more directly, however, 
higher-resolution observations at submillimeter wavelengths are essential. 

This paper presents simultaneous observations of dust continuum emission at 340 GHz and 
CO (J = 3 — 2) from FN Tau with the Submillimeter Array (SMAfl. High angular resolution 
provided by the SMA has allowed us to constrain better the distribution of circumstellar material. 
We describe the details of the observations in §2 and present the results in §3. Based on our results 
as well as those obtained by previous studies, we discuss in §4 the circumstellar structure of FN 
Tau, including the disk parameters and the origin of the near-infrared nebulosity. 



Observations 



SMA observations of FN Tau were carried out on 2008 July 11 (UT). 4K-cooled SIS mixer 
receivers were tuned at the central frequencies 336.00488 GHz and 345.53059 GHz for the lower 
and upper sidebands, respectively, and the continuum emission was received with a bandwidth of 
1.7 GHz for each sideband. Data of CO ( J = 3 — 2) were simultaneously obtained with a frequency 
resolution of 101.5 kHz, corresponding to a velocity resolution of 0.088 km s^^. The extended 
antenna configuration consisting of 7 antennas provided baselines ranging from 39 kA to 258 kA. 
The largest brightness distribution to which these observations a re sensitive enough (rn ore than 
50% level) is ~ 2.3" in FWHM when the distribution is Gaussian (jWilner fc Welchlll994l ). Uranus 
was observed to obtain the passband characteristics of the system, while the two quasars 3C111 
and 0336-1-323 were observed every 15 minutes to track the time variation of the complex gain. The 
system temperature during the observations, Tgys, was 230 — 450 K. 

Visibility data were calibrated with the IDL MIR package, which was originally deve lop ed 
by the Owens Valley Radio Observatory as the MMA software package (jScoville et al.lll993l ) and 
later adapted for the SMA. After the correction for the visibility amplitude using Tgyg, the gain 
table constructed from the 3C111 data was applied to the visibilities of FN Tau and 0336-1-323; 



^ The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the Academia 
Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and the Academia Sinica. 



- 5 - 



the flux density of 3C111 at 340 GHz was derived to be 1.04 Jy from the comparison with Uranus. 
Calibrated visibilities of 0336+323 were carefully inspected to identify visibilities taken under poor 
conditions. All the visibilities identified were flagged out before further analyses in §3. Data 
from both the sidebands were combined in the analyses of the continuum emission, hence the 
effective central frequency is 340.767735 GHz. Imaging from the calibrated visibilities was made 
with Astronomical Image Processing System (AIPS) developed at the National Radio Astronomy 
Observatory (NRAO). 



3. Results 

3.1. Continuum Emission 

Figure [1] shows the maps of the continuum emission from FN Tau constructed with the CLEAN 
algorithm under two different visibility weightings. Natural weighting provides a synthesized beam 
of 0'.'78 X 0'.'71 (P.A. = 38°) in FWHM, and the rms level in the emission-free regions is 1.7 mJy 
beam~^, corresponding to 0.0324 K in brightness temperature. Gaussian fit to the detected emission 
reveals that the emission has a total flux density of 36.7 it 2.9 mJy with no significant extension 
compared to the beam size. A sharper beam, 0'.'74 x 0'.'57 (P.A. = 68°) in FWHM, is given in the 
uniform weighting, but the detected emission is again point-like with a flux density of 31.5 it 4.7 
mJy. Although these results can be regarded as the detection of a point source, slightly smaller 
flux density with the uniform weighting may suggest that the spatial extent of the emission could 
marginally be resolved. Similar comparison has also been made for 0336+323, which must be a 
perfect point source, and the detected flux densities are 235.3 ± 8.2 mJy for the natural weighting 
and 233.8 ± 7.8 mJy for the uniform weighting; unlike the case of FN Tau, there is little difference 
between these two weightings. 

These situations can also be confirmed by the visibility plots (Figure [2]) . When a point source 
located at the phase center is observed, the real part of its visibility is a constant equal to its flux 
density while the imaginary part is zero for all the baselines. Indeed, these characteristics can 
be clearly seen in the visibility plots for 0336+323. Visibility plots for FN Tau also show similar 
characteristics, but the imaginary part at < 165 kA is nearly zero while the real part in (135—165) kA 
is slightly smaller than that at shorter baselines. These can be explained by a circularly symmetric 
distribution of the emission with a finite extent. The submillimeter continuum emission from FN 
Tau is likely to have such a nature because t he circumstellar ne bulosity at near-infrared exhibits a 



circularly symmetric brightness distribution (iKudo et al.l 120081 ). At baseline lengths greater than 



180 kA, however, the imaginary part has non-zero values, possibly due to higher phase noise. Thus 
the deviation from the perfect point-source case cannot confidently be identified in our results. 



Upper limit to the apparent size of the emission is estimated from a comparison with the 
analytic solution for the case of uniform brightness. Visibility amplitude for a disk with uniform 
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brightness, |1^(Z)a)|, is expressed by 



ttDx9 



(1) 



where Ji is Bessel function of the first kind, Fq and 9 are the total flux density and apparent radius 
of th e disk, respectively, an d Dx is projected baseline length in units of the observed wavelength 
(e.g., 



Wilner k Welchlll994l ). \V{Dx) \ peaks at Dx 



with Fq and takes the 80% level of Fq when 
0.21. In the case of FN Tau, the real part of the visibility is ~ 30 mJy, or 80% of the 
total flux (37 mJy), at Dx ~ 150 kA. This means 6 = 0.29", or 41 AU in radius, when Equation 
([1]) is applied. Indeed, the solution for 6 = 0.29" still traces the SMA data fairly well as shown 
in Figure [2^a). The outer radius of the detected emission is therefore estimated to be < 41 AU, 
significantly sma ller than that of the near- infrared nebulosity (r = 260 AU for the bright part; 
Kudo et al.ll2008l ;i. 



Continuum flux density of FN Tau was also measured to be 31 it 9 mJy at 230 GHz with the 
11" beam of the IRAM 30- meter telescope (jBeckwith et al.lll99d ). Flux density measured at 340 
GHz with the SMA (37 mJy) is significantly smaller than the value simply extrapolated from the 
measurement at 230 GHz with F,^ oc i/^"^^ , where /3 is th e exponent of dust ab s orption coefficient for 



the d isks around T Tauri stars ranging from to 2 (e.g., iKitamura et al 



2002 



Andrews &: Williams 



20071 ). This suggests that some fraction of dust emission is missed in the SMA observations; this 
issue will be examined further in M.ll 



3.2. CO Emission 



After averaging visibilities in five contiguous frequency channels, we made maps of CO (J = 
3 — 2) with a velocity resolution of 0.44 km s~^. Natural weighting provides an rms noise level 
of 100 mJy beam~^ with a synthesized beam of 0'.'78 x 0'.'71 (P.A. = 38°), or 2 K in brightness 
temperature, but no significant emission was detected. 

The disk component which is < 41 AU in radius and is detected with the dust continuum can 
also emit the radiation of CO (J = 3 — 2), but our sensitivity seems insufficient to reveal the CO 
emission as follows. When the disk is perfectly face-on and emits the CO line whose brightness 
temperature is ~ 19 K, as suggested by the model fit in ^4.11 beam dilution will make the observed 
intensity as low as 7.7 K, which is just below the 4o" level. Note that this estimate is for an 
optimistic case because the emission from the di sk is assumed to fall in a single velocity channel 
because of the perfect face-on configuration (e.g.. lOmodaka et al.lll992l ). The Keplarian rotational 
velocity around FN Tau will be ~ 1.5 km s~^ at r = 41 A U, and the observ ed velocity width will 
be > 1 km s~^ if the inclination angle is 20°, as inferred by lKudo et al.l ( 20081 ). In such a case, only 
a partial area will emit the CO radiation for a single velocity channel and the observed intensity 
will become lower. Further niore, FN Tau is loc ated in the L1495 cloud where the extended CO 
emission is very strong (e.g., iMizuno et al.lll995l ) while the SMA is sensitive only to the increment 
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from these extended components, which makes the detection even more difficult. More sensitive 
observations will be required to reveal the compact CO emission associated with the star. 



4. Discussion 



Tau 



4.1. Fitting of Compact Dust Emission with a Power-law Disk Model 

4-. 1.1. /3 = 1 case 

T he continuum emis sion detected with the SMA should originate from the disk around FN 
pOOSi ) adopted a power-law disk model with inner and outer cutoffs (e.g., 



Kudo et al 



Kitamura et al 



2OO2I ) and derived disk parameters by fittin g the spectral energy distribution, in- 



cluding measurement at 230 GHz with the IRAM telescope (jBeckwith et al.lll99Cll ). They assumed 
that the gas-to-dust mass ratio is 100 and that the mass opacity of dust particles, k,/, at millimeter 
and submillimeter wavelengths can be expressed as 



0.1 



1012 Hz 



cm g 



(2) 



with /3 = 1. iKudo et al.l (|2008l ) also assumed that the bright inner part of the near- infrared 



nebulosity, 3.7" in diameter, coincides with the disk size. The emission component detected with 
the SMA at 340 GHz, h o wever , is much weaker than 100 mJy expected from the disk parameters 
obtained by iKudo et al.l ([200a). Moreover, its size is so small that the ~ 0.7" beam of the SMA 
could not resolve the spatial distribution. This is not due to an artifact c aused by the finite 
sensitivity of our observations; as shown in Figure O the model disk derived by Kudo et al. (j2008l ) 
has a brightness temperature at 340 GHz higher than the 3<t level of our observations inside ~ 100 
AU in radius, and the SMA would resolve the emission if the power-law disk h ad a radius greate r 
than 100 AU. The present observations, therefore, prove that the assumption bv lKudo et al.l (|2008l ) 
that the radius of the power-law disk coincides with the extent of near-infrared nebulosity is invalid. 

We therefore reexamine the disk parameters by fittin g the spectral energy distribution with the 
same power-law disk model as that in lKudo et al. (2008) but under a different assumption that the 
outer radius of the disk is 41 AU, or the upper limit to the outer radius of the emission detected 
with the SMA (se e §3. ID . Indeed, the emission at A < 60 /^m expected from the model disk in 
Kudo et al.l (120081 ) is attributed to the region well inside 41 AU in radius. In the present model fit, 
we include the flux d ensity at 340 GHz wi th the SMA, but the result at 230 GHz obtained with 
the IRAM telescope (jBeckwith et al.lll99d ) is excluded because its beam size (11" in FWHM) is 
so large that emission from a cold extended component outside the power-law disk may also be 
contained. At infrared wavelengths, photometric data obtained with the Spitzer Space telescope 
at 3.6, 4.5, 5.8, 8.0, 24 and 70 /im (Spitzer Taurus Catalog October 2008 v2.ljl are included in 



^The photometric data in the Spitzer Taurus Catalog are taken from 'http://irsa.ipac.caltech.edu/index.htmH 



- 8 - 



addition to those used by iKudo et al.l (j2008l ). while the data point at 100 //m is excluded because 
this is merely an upper limit obtained with IRAS. The best fit disk parameters are summarized in 
the first row of Table [H and the comparison between the spectral energy distribution of the model 
and the ob served phot o metri c data is shown in Figure \M,a). Compared to the disk parameters 
inferred by iKudo et al.l (| 20081 ) , the surface density distribution in Table [1] is comparable while the 
temperature is slightly lower with a steeper radial distribution. Bec ause of the s i naller outer radius, 
the total disk mass is 1.1 x 10~'^ Mq, only 15% of the estimate by lKudo et al.l (j2008l ). 



Although the result shown in Figure (H^a) is for the best fit case, it still exhibits notable 
discrepancies in two wavelength ranges: (i) at near-infrared of A = 0.8 — 2.2 ^um, in which the 
photometric data are mainly taken from the 2MASS catalogue, and (ii) at wavelengths longer 
than 5.8 /xm. As shown in the following, the discrepancy at near- infrared is probably due to 
the contribution of the circumstellar nebulosity in scattering light, which is not included in the 
model. Photometric magnitude of FN Tau in the H-band from the 2MASS catalogue is 8.67 mag, 
corresponding to 6.4 x 10~^^ erg cm~^ in XFx, but this contains both the thermal radiation from 
the star-disk system and the scattering light from the circumstellar nebulosity. The contribution 
from the star in the model is calculated to be 3.7 x 1 0"^*^ erg s~^ cm~^, which is consistent with 
photometric measurements of M5 stars (|Pickleslll998l ). Thermal radiation from the model disk is 
estimated to be ~ 10% of the stellar radiation at this wavelength, hence the amount of discrepancy, 
or the contribution of scattering light, is derived to be ~ 2.3 x l O"^'^ erg s~ ^ cm~^. On the other 



hand, the flux density of the nebulosity in the H-band image bv lKudo et al.l (|2008l ) is 6.9 mJy, or 
1.3 X 10^^*^ erg s~^ cm~^. This is only ~ 55 % of the excess emission estimated above, but it does not 
include the contribution from the area blocked by the occulting mask of the coronagraph (~ 0.8" in 
diameter) and the spider pattern of the telescope. If the contribution from the blocked area is taken 
into consideration, the nebulosity component can account for the discrepancy at near-infrared in 
the model fit. 

The situation at longer wavelengths, which are important in deriving the disk temperature dis- 
tribution, seems more complicated. Figure [D^a) shows that the observed fluxes in A = 5.8 — 10 /im 
are systematically smaller than the model while those in A = 12 — 60 fim are systematically larger 
than the model. These disagreements could be resolved if we adopt a more elaborat e disk model 
with a flared surface layer of superheated dust grains (e.g.. IChiang &: Goldreighlll997l ) because the 
surface layer can contribute additional fluxes mainly at the latter wavelength range. We should 
examine in detail these effects when we would investigate the vertical structure at smaller radii, 
but this is beyond the scope of this paper. To estimate the total mass or the surface density of the 
disk based on our SMA results, a simple power-law disk model may be acceptable because almost 
whole the disk emission at millimeter and submillimeter wavelengths are expected to originate 
from a single layer, i.e., the interior region near the mid-plane which is the dominant mass reservoir 
(jChiang Sz Goldreighl 119971 ). For the temperature distribution of the interior region, a better esti- 
mate could be obtained with a model fit in which photometric data that might be severely affected 
by the contribution from the surface layer are eliminated. As shown in Figure \S(b) and the fourth 
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row of Table dl the result of the model fit without data at A = 12 — 60 /um, which exhibit significant 
upward deviations from the model in Figure HH^a) , gives us a disk mass 40 % higher than the case 
shown in Figure [l^a) due to the lower derived temperature. For the /3 = and /3 = 2 cases in 
^4.1.21 and in Table [H we will only show the results of the model fits in which all the infrared data 
are included, but the disk mass can be higher than those shown in Table [1] by 40 % in the model 
fits without data at A = 12 — 60 ^m, as similar to the (3 = 1 case. 

The disk models in Figures IHa) and \Mb) predict that the flux density at 230 GHz is 13 it 3 
mJy. This is significan tly smaller than the flux density measured by the 11" beam of the IRAM 



telescope, 31 it 9 mJy (jBeckwith et al.lll990l ). Since this measurement was made with a chopping 
position for the sky subtraction just 30" (4200 AU) apart from the star, contributions by the 
extended background should mostly be canceled out and the detected emission is likely to originate 
from some circumstellar components. One can also confirm easily that the contamination from the 
extended CO ( J = 2 — 1) line is negligible when one takes into account the fact that the frequency 
bandwidth of this observation was 50 GHz. The above comparison, therefore, suggests that the 
fiux density measured with the IRAM telescope contains a slight excess of 18 it 9 mJy from the 
power-law disk model with /3 = 1, and that there is another component which is missed in the 
present SMA observations. The possible origin of this missing component, together with apparent 
disagreement in size between the emission seen in the SMA image and the near-infrared nebulosity, 
will be discussed in N4.2[ 



4- 1-2. Cases of Other (3 Values 



Since there is no photometric data at other millimeter and submillimeter wavelengths with a 
similar beam size as that of the SMA, we cannot directly determine j3 inside the disk, or r = 41 
AU. We therefore made model fits with other (3 values, and 2, to evaluate uncertainties in derived 
physical quantities. The obtained disk parameters are summarized in Table [U and the comparisons 
between the spectral energy distributions of the models and the observed photometric data are 
shown in Figure [5l The parameters related to the temperature distribution in Table [1] are almost 
insensitive to the change in (3 because these are essentially determined by data points in the infrared 
range where the disk is optically thick. The surface density and disk mass, on the other hand, can 
be varied by a factor of 3 compared to the (3 = 1 case according as the change of K,y following 
Equation ([2]). 

Although there are large uncertainties in the estimates for the disk mass and surface density, we 
can obtain several robust conclusions. First, the estimated disk mass is at most 5.9x 10~^ -^©j which 
corresponds to 1.4 times the maximum mass in the (3 = 2 case of Table [H Thi s is .just an order o f 
Jupiter mass and seems insufficient to build a gas giant planet. As discussed bvlKudo et al.l (|2008l ). 
it is most likely that only small terrestrial or icy planets will form (jKokubo &: Idall2002l ). Secondly, 
the disk size (r < 41 AU) is significantly smaller than the typical disk radius for solar-type T Tauri 



stars derived from a similar model fit including SMA data (r ~ 200 AU; lAndrews &: Williams 
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20071 ). but the spectral energy distribution of FN Tau can be explained by a power-law disk model 



as sim ilar to the cases of other solar-type T Tauri stars (jKitamura et al.ll2002l : [Andrews &: Williams 
20071 ). Together with the fact that the disk ma ss relative to the star (0.3 — 3 %) is within the range 



of the statistics for solar- type T Tauri stars ([Andrews &: Williamsl 120071 ). FN Tau has probably 
been formed through the standard process of star formation, i.e., the gravitational contraction 
of a cloud core, but not through dyn amical ejection from a gravitationally unstable circumstellar 
disk (jBate. Bonnell. &: Bromm |20oi). Thirdly, all the best-fit model disks whose parameters are 
listed in T able [J cannot fully re produce the flux density at A = 1.3 mm obtained with the IRAM 
telescope (jBeckwith et al.lll99d ). indicating that the two results with the SMA and the IRAM 
telescope cannot simultaneously be explained by a power-law disk only. This implies that there is 
another more extended component that cannot be well described by a power-law disk model and 
is missed in the SMA observations. 



4.2. Nature of the Missing Extended Component 

As shown in §4.1^ the best-fit power-law disk model cannot fully reproduce the flux density at 
230 GHz obtained with the IRAM telescope in all the cases of /3 = — 2. This implies that, in 
addition to the compact power-law component, there is a more extended component that is missed 
in the SMA observations; we will denote this component by "halo" in the following. 

We can roughly evaluate the lower limit to the size of the halo by requiring it not to be detected 
with our SMA observations. In the case of /3 = 1, for example, the flux density of the model disk 
component at 230 GHz is 13 mJy (see Tabled]), hence the median for the expected contribution of 
the halo at this frequency is 18 mJy. Assuming that dust particles in the halo have the same /3 as 
those in the disk, the flux density of the halo component at 340 GHz is estimated to be 58 mJy. If 
the emission of the halo is uniformly distributed inside 2.61" (i.e., 366 AU) in radius, it will have 
a brightness temperature equal to the la level of the map in Figure [11(a) and cannot be detected. 
This estimate can also be verified in the visibility plot; the contribution of the halo at Dx = 45 
kA will be 3.5 mJy (i.e., 6% of the total flux density of the halo) according to Equation ([T|) with 
9 = 2.61", which is much smaller than the standard error of the visibility amplitude (see Figure [2]). 
We have summarized in Table [2] the lower limits to the size of the halo for other /3 cases including 
their uncertainties. 

The estimated lower limits to the halo's radius in Table [2] are between 174 AU and 574 AU, 
depending on the assumed /3 value. All of these numbers, however, are smaller than the beam size 
of the IRAM telescope, thus its photometric measurement can completely contain all the emission 
from the halo. More importantly, the estimated size of the halo is quite reminiscent of the size of 
the near-infrared nebulosity: 1.85" (260 AU) in radius for the brighter part where the brightness 
distribution S{r) f ollows S(r) pa r~ ^'^ and 2.85" (400 AU) in radius for the tenuous periphery 



with S{r) oc r'''-^ (|Kudo et al.ll2008l V These comparisons suggest that the circumstellar structure 



around FN Tau may consist of two components: a compact power-law disk (r < 41 AU) and a 
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surrounding halo that probably corresponds to the origin of the near-infrared nebulosity. The SMA 
observations have revealed only the former component while both the components were detected 
by the IRAM telescope. 

When all the infrared data are used in a model fit, the best-fit parameters of the disk give 
T 19 K at the outer edge of the disk, rout = 41 AU, regardless of the assumed /3 (see Table [I]). If 
we assume that the temperature in the halo, Thaio; is also the same as that at rout and uniformly 19 
K, we can estimate its mass (Mhaio) and average surface density (Shaio) by the following equations: 



3.8 X W'^Mq 



d 



f -^340 GHz 

V 58 mJy ) \im pc ) VO-0341 cm^ g-i 



^halo 



-^halo n 1 n-2 -2 

— 2 — = 7.9 X 10 g cm 

halo 



vrr 



My 



halo 



?^halo 

3.8 X IO-^Mq/ V366 AU 



-1 



-2 



(3) 
(4) 



Mhaio and Shaio for all the cases derived by Equation (l3|) or (jl]) are shown in Table [21 in which Shaio 
are normalized by the disk surface density at rout (see Table [T|). Note that Mhaio in Table [2] will be 
higher by ~ 40 % when we adopt the disk parameters obtained by a model fit without photometric 
data at A = 12 — 60 ^m (see ^4.1.ip . but the mass and surface density of the halo relative to the 
inner disk remain unchanged because the disk mass and its surface density will also be higher by 
the same factor. Table [2] clearly indicates that the halo can contain mass comparable to or at most 
8 times greater than that of the inner power-law disk, but Shaio should be lower than the surface 
density of the inner disk by an order of magnitude. 

Based on the fact that the surf ace brightne s s distr ibution of the near-infrared nebulosity, S{r), 



follows the relation S{r) oc r~^'^, iKudo et al.l ( 20081 ) concluded that its origin is the scattering 
light at the surface of a flared disk. The new insight obtained with the present SMA observa- 
tions, however, is that the circumstellar component outside at least 41 AU in radius cannot be 
described well in terms of a power-law disk model. Indeed, a power-law fo rm of S{r) can also 



be explained either by an optically thin, spherically symmetric envelope (e.g., iTamura et al. 



1988 



Weintraub et al. 



19921') or by a flattened envelope with larger optica l depth in almost pole-on con- 



figuration (jWhitney &: Hartmannlll993l ). though iKudo et al.l (j2008l ) did not examine these possi- 
bilities in detail. Among these possibilities, an optically thin and spherical envelope is unlikely; 
Ehaio = 7.9 X 10^^ g cm~^ in Equation (HI) would correspond to Ay ^ 9. 4 mag to the star in the 
usual relation for the interstellar medium (jBohlin. Savage. &: Drake! Il978l ) when a spherical enve- 



lope is assumed, but this is much larger than the actual Ay to the star (1.35 mag; iMuzeroUe et al 
2OO3I ). 



What is then the nature of the halo ? One possibility is that this is part of a flared circumstellar 
disk but the periphery regions that cannot be described well by a simple power-law disk model. 
Recently, iHughes et al.l (j2008l ) investigated the apparent discrepancy between gas and dust outer 
radii derived from millimeter and submillimeter observations of protoplanetary disks. They found 
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that a disk model described by power laws in surface density and temperature that are truncated 
at an outer radius does not simultaneously reproduce the gas and dust emissions, but a simple 
accretion disk model that includes a tapered exponential edge in the surface density distribution 
does well. This is because the surface density in the outer tapered regions can be so low that the dust 
continuum emission at submillimeter is negligible but large enough to produce detectable emission 
in a probe of higher emission coefficient such as a rotationa l traii sition of CO or scattering light 
at near-infrared. In fact, one of the sample in lHughes et al.l ()2008l ). TW Hya, shows quite similar 
features to those seen in FN Tau; it exhibits near-infrared nebulosity whose size (r = 170 AU) is 
significantly larger than the size of the submillimeter continuum emission determined by visibility 
fitting with a power-law disk model (r = 60 AU), and the radial brightness d istribution of the 



nebulosity has a power-law form which is almost the same as that around FN Tau (jWeinberger et al 



2OO2I ). The inferred surface density of the halo (see Equation @ and Table [2]), which is significantly 



lower than the inner power-law disk, is qualitatively similar to the outer regions of the model disk 
with an exponential tail. Another possibility is that the halo is a diffuse remnant of the protostellar 
envelope wit h flattened structure (Whitney &: Hartmannlll993l ). The age of FN Tau is estimated to 
be < 10^ yr ( Muzerolle et al.ll2003l ). hence it is not surprising that FN Tau is still accompanied by 
such a component. Although these two candidates for the halo or the origin of the nebulosity are 
quite different in nature, it is unclear which is the case because both of them will have a scattering 
surface that can account for the radial brightness profile of the near-infrared nebulosity. Further 
sensitive observations in line emission will be critical to understand the nature of the halo. 
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Fig. 1.— Continuum emission at i/ = 340.768 GHz toward FN Tau with the SMA. (a) Map with 
natural weighting. The synthesized beam is (/.'TSxC/.'?! and P. A. = 38° in FWHM, and the rms noise 
level is 1.7 mJy beam~^. (b) Map with uniform weighting. The synthesized beam is (y.'74 x 0f.'57 
and P. A. = 68° in FWHM, and the rms noise level is 2.9 mJy beam~^. Contour interval is 2a, 
starting at ±20" levels. Dashed lines indicate negative levels. Ellipse at the bottom left corner of 
each panel shows the beam size in FWHM. Both axes are measured from the phase center (a, 5) = 
(04^14™14?59, 28°27'58") in J2000. 
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Fig. 2. — Visibility plots for FN Tau (a, b) and 0336-1-323 (c, d). Both the real and imaginary parts 
of the visibilities are binned and averaged every 30 kA. Standard errors are also plotted. For FN 
Tau, the positional offset of the emisson peak from the phase center, (—0.2", -1-0.2"), is corrected. 
No correction was made for 0336-1-323 because the emission peak coincides with the phase center. 
Dotted line in (a) indicates the analytic solution for the disk with uniform brightness whose total 
flux density is 37 mJy and the apparent radius is 0.29". 
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Fig. 3. — Brightness distribution at u = 340 GHz of the model disk in iKudo et alJ (|2008l ). Dotted 
hne indicates the 3a level of the continuum map with natural weighting, 0.097K (see ^3.1|) . 
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Fig. 4. — Spectral energy distribution of FN Tau and that of a model disk which is best-fitted to 
(a) all the infrared data and (6) infrared data without A = (12 — 60) fxm under the assumption of 
/? = 1 . The parameters of both the model disks are listed in Table [H Gray dashed lines indicate 
the contribution of the stellar radiation. The data points indicated by gray squares are not used 
in the model fit. Error bars for data at A = 70 fim (20 %) and at 1.3 mm (30 %) are shown, but 
those for other data are not shown because each error is so small (< 10 %) or is unknown. 
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Fig. 5. — Model fit of tlie spectral energy distribution of FN Tau by a power-law disk model with 
different /3 values; the (3 = 1 case is shown by the solid line, the /3 = case by the dotted line, 
and the (3 = 2 case by the dashed line. The parameters of all the model disks are listed in Table 
[TJ Gray dashed line indicates the contribution of the stellar radiation. The data point at A = 1.3 
mm, indicated by a gray square, is not used in the model fit. 
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Table 1: Best Fit Parameters for a Model Disk 







TiAV 


Q 








-^230GHz 






(K) 




(g cm~^) 


(^0) 




(mJy) 


/3 = 


lb 


141 ±3 


-0.543 ± 0.011 


38 ±9 


(1.11 ±0.27) X 10- 


-3 


13 ±3 


/3 = 




144 ±3 


-0.534 ±0.012 


11±3 


(3.23 ± 0.87) X 10" 


-4 


16 ±3 


P = 


2^ 


141 ±3 


-0.545 ±0.011 


115 ±28 


(3.37 ±0.83) X 10" 


-3 


9±2 


P = 


1^ 


126 ±5 


-0.577 ±0.015 


53± 14 


(1.55 ± 0.40) X 10" 


-3 


13 ±3 



"Gas-to-dust mass ratio is assumed to be 100. 

''Fitting results with all the infrared data. 

■^Fitting results without using data at A = (12 — 60) /xm. 

Note. — All the uncertainties in derived physical quantities are evaluated so that each photometric measurement 
contains at most 20% uncertainty, which is the case for the data point at A = 70/^m; Actual uncertainties for other 
data points at other wavelengths are much smaller than 20%. Fixed parameters are as follows: i = 20°, rin = 0.015AU, 
rout = 41AU, r,tar = 0.0105AU, p = 1, T.tar = 3240K, and Ay = 1.35mag. 
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Tahle 2: EsUmalcd Quanlilics for the Halo 





-^230GHz 
(mJy) 


-^^340GHz 
(mJy) 


Minimum Radius 

(") 


Mhalo^ 

(Mo) 


Mhalo/Mdisk 


^halo/^disk'^ 


p = 


1 


18 ±9 


58 ±29 


2.6ii°:i 


(3.8 ± 1.9) X 10-3 


3.4 ± 1.7 


< 0.086 


p = 





15 lb 9 


33 ±20 


i.96i°:i 


(7.3 ±4.4) X 10"^ 


2.2 ±1.3 


< 0.10 


p = 


2 


21 ±9 


100 ± 43 


3.43l°|I 


(1.9 ±0.8) X 10-2 


5.7 ±2.4 


< 0.083 



"Assuming that the temperature of the halo to be 19 K and that the gas-to-dust mass ratio to be 100. 
''Surface density at the outer radius (= 41 AU) is adopted as Ti^iak- 



